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Abstract 



^3 ' The thermodynamical partition function of the Duffm-Kemmer-Petiau theory is 



evaluated using the imaginary-time formalism of quantum field theory at finite tem- 
perature and path integral methods. The DKP partition function displays two features: 
(i) full equivalence with the partition function for charged scalar particles and charged 
\ massive spin 1 particles; and (ii) the zero mode sector which is essential to reproduce 

^ ■ the well-known relativistic Bose-Einstein condensation for both theories. 
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1 Introduction 



Bose-Einstein condensation ^ |2] (BEC) is a phenomenon long associated to the liquid 
helium studies ( 4 He and 3 He- 4 He mixtures). However, recent research deals with BEC 
| shown that it touches several areas of modern physics 3 : in thermodynamics j2| BEC 

occurs as a phase transition from gas to a new state of matter, quantum mechanics view 
BEC as a matter-wave coherence jSJ arising from the overlapping de Broglie waves of the 
atoms and draw an analogy between conventional and "atom lasers" fBj , quantum statistics 
explain BEC as more than one atom sharing a phase space cell, in the quantum theory of 
atomic traps [7] many atoms condense to the ground state of the trap, in quantum field 
theory BEC is commonly related to spontaneous symmetry breaking UJ EI • 

In quantum field theory BEC the complex scalar field was and is used for studying the 
thermodynamical properties of physical systems composite of bosonic particles with spin 0. 
However, at zero temperature, there is an alternative way to study not only the properties 
for complex scalar field but also the complex vectorial field that it is known as the massive 
Duffm-Kemmer-Petiau (DKP) theory [H] . 

At zero temperature, one important question concerning DKP theory is about the equiv- 
alence or not between its spin and 1 sectors and the theories based on the second order 
Klein-Gordon (KG) and Proca equations, respectively. From the beginning of the 50's the 
belief on this equivalence was perhaps the principal reason for the abandon of DKP equation 
in favor of KGF and Proca ones. However, in the 70's this supposed equivalence began to be 
investigated in several situations involving breaking of symmetries and hadronic processes, 
showing that in some cases DKP and KG theories can give different results (for a historical 
review of the development of DKP theory until the decade of 70's see the reference \V2\). 
Moreover, DKP theory follows to be richer than the KG one with respect to the introduction 
of interactions. In this context, alternative DKP-based models were proposed for the study 
of meson-nucleus interactions, yielding a better adjustment to the experimental data when 
compared to the KG-based theory JSj- In the same direction, and guided for the formal 
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analogy with Dirac equation, approximation techniques formerly developed in the context 
of nucleon-nucleus scattering were generalized, giving a good description for experimental 
data of meson-nucleus scattering. The deuteron-nucleus scattering was also studied using 
DKP, motivated by the fact that this theory suggests a spin 1 structure from combining two 
spin | [H[TS]. 

Recently there have been a renewed interest in DKP theory. For instance, it has been 
studied in the context of QCD j!6) . covariant hamiltonian dynamics ^J], in the causal 
approach |18j . in the context of five-dimensional galilean covariance |19| . in the scattering 
if+-nucleus [201, in curved space-times |23 ED > etc. 

These examples, among others in literature, in some cases break the equivalence between 
the theories based on DKP and KG/Proca equations, such as in or in Riemann-Cartan 
space-times [23 El- Nevertheless, the question about the equivalence or not still lacks a 
complete answer nowadays. 

As above mentioned, all accomplished studies on the DKP theory were made at zero 
temperature. The aim of this work is to study the thermodynamics and Bose-Einstein 
condensation for the massive charged particles by using the DKP theory. In section 2 we 
present the massive DKP theory in the Minkowski space-time and we make a resume of the 
constraint analysis of the model which can be seen in |25j . where it has been shown that 
all the photon Green's functions coincide in DKP and KGF theories. In the section 3 we 
introduce the partition function and the generating functional of the correlation functions 
of the theory, and we specialize, separately, the spin and spin 1 sector for the explicit 
calculation of the respective partition functions. We also analyze explicitly the zero mode 
contribution in both cases. In the section 4 we gives our conclusions and perspectives. 



2 The Duffin-Kemmer-Petiau theory 

The Duffin-Kemmer-Petiau equation in Minkowski space-time is given by 

i^d^ %p - m%p = 0, (1) 
where the matrices (3^ obey the DKP algebra, 

P*pPfi> + /3 u (3"/3 a = p*if + ^V" (2) 

with if 1 " being the metric tensor of Minkowski space-time with signature (H ). The 

/3 M are singular matrices which have only three irreducible representations of dimensions 1, 
5 and 10. The first one is trivial, having no physical meaning and the other two correspond 
to fields of spin and 1, respectively. 

The DKP equation given in is obtained from the following Lagrangian density 

C = I i>^d^ - l - drfpt - , (3) 

where if) = ip^rf and rf = 2(/3 ) 2 — 1, with {if) 2 = 1 

In , the Hamiltonian canonical approach to the DKP theory is developed in the 
components and the matrix forms. We make a resume of the matrix form of this analysis 
to follow. The canonical conjugate momenta are 

The Hamiltonian density is given by 

W = - \ Wdki) + \ d^{i k i> + m$i> . (5) 
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There is a set of second-class constraints 



(6) 



cj = M(t /3 fc a fc ^ - mVO , w = (t d k ip{3 k + mlp) M , (7) 
where M = 1 — (f3 ) 2 is a singular matrix. 

3 The Partition function 

The massive DKP Lagrangian density © has a global U(l) symmetry 

^ -> = e ia il) 

i> -> i 1 = e- ioi $ (8) 
whose conserved current is j M = ipP^tpi an d the respective conserved charge 

Q = Jd 3 xf(x) . (9) 

Then the partition function in the imaginary-time formalism for the massive DKP field is 



Z = N(J3) I V P ,p Vp-- J Vip DV 5 (9) S (9) S (u) S (sj) 

periodic 



(10) 



exp < dr d x 



ipj, d T ip + id T Tpp^- H + fj,j° 



where N(f3) is an infinite normalizing factor which will be determined later 26 . We integrate 
over periodic DKP field due to its bosonic character. 

To compute the correlation functions of the theory we couple the DKP field to external 
sources fj and r], thus the generating functional is read 



Z[rj, ff] = N(P) I Vp^ Vp^ I V^V^S (9) 6 (9) 6 (w) S (uJ) 

periodic 



(11) 



exp \ dr d x 



ip^ d T tp + id T tp p^ - H + fij +r)ij) + tpr) 



After the integration on the canonical conjugate momenta, we define the matrix (3 T = i/3 ; 
and the constraints to and uJ can be rewritten as 

u = MDip , ZJ = ThpM 

due the property /3°M — = M/3°, and we have used a short notation, 

D?/> = a d a tp -mtp , r>0 = idjpp" + imp, (12) 

where a — r, 1, 2, 3. With all definitions noted above, the generating functional appears 
to be 



Z[r},rj] = iV(/3) J VijiVijj S(MBrp)S (D^M) e f 



(13) 



periodic 
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where the action S is given by 



S = I dr I d 3 x ( ~ tp Dtp - i Dip ip + fj,ipp3°ip + r]tp + iprj 



(14) 



After the exponentiation of the <5-functions, we use the periodicity condition over the fields 
to make some integration by parts, and get the following expression for the generating 
functional 



Z[r),Tj] =N(J3) I ViPViP VpVW e sl 

periodic 



(15) 



where the auxiliary fields p and a were used to exponentiate the constraints; and the action 
S 1 is given by 



S 1 = / dr d 



(16) 



(17) 



ip (D + p(3") ip - i DMcr ip-iip DMp + rjip + ^prj 
Into the generating functional l|15|) . we make the following change of variables, 

Ip — ► Ip + 1p s , tp —* 4> + i>s 

such the fields ip s and ip s satisfy the following equations 

(D+z/ji-^Mp + i) = ; -Blp~ s + plp s (3 -iDM^+rj = 0. 
And after it, we carry out another change of variables on the auxiliary fields p and c, 

a^a + a s , p^p + p s (18) 
where the fields p s and a s satisfy the equations 

mMp s - iM-q = ; mW s M + irjM = . 
Thus, after some additional calculations in the generating functional in l|15f) . we get 



Z[r),rj] =cxpj-^ drjd 3 x fj (D + pP°) \ + i ffMr]^ Z, 



(19) 



where Z is the partition function given in l|lUfl which, after all operations made above, it is 
read as 



Z = N(/3) / VipVip VpVa exp 
criodic 



per 



J drj rf 3 x ip(D + npf)tp + moMp 



(20) 



We can note that the functional integration over the fields ef and p give an infinity quan- 
tity, due to singular character of the M-matrix, which is irrelevant to the thermodynamical 
properties of the physical system. Thus, we can write the partition function as 



Z = N'N(fJ) [ VipViP exp 

periodic 



f drJd 3 xtP(D + pf3°)ijj 



(21) 



The fields may be expanded as a Fourier series (2HJ 

d 3 p 

(2tt)3 



1p(T,X) = Ipzero+Jj^j- 



1pn(p) 



(22) 



ip(r,: 



d 3 p 



7 iv/i] 

Vzero + p 2^ J ( 27r 



e i(p~+«„rty n(p) 
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where ipzero and ip zero are independent of (r, x) and carry the full infrared character of the 
DKP field; that is, VVi=o(p = 0) = and ip n=0 (p = 0) = 0. And due to bosonic character 
of DKP field, we impose the constraint of periodicity such ip(r = 0, x) = ip(r = (3, x) for all 
x, thus the frequency will be uj n = 2-7m//3, n = ±f , ±2, .... 
Substituting (J2U) into (|2TT) . we find 

p r _ 

■ ■ / d 3 x V(D + fj./3°)ip = f3Vi/j zero {-m + M/3°)VWo + 5' , (23) 



where the free infrared action is 

S' = ~-^J^^- n {-Y>)^P°^n + i^) + P k p k + m }Mv) ■ (24) 

The partition function given in (|2f |) . can be rewritten using the equation l|24[l as 

Z = N'N(p) Z zero Z' , (25) 
where Z zero is the zero mode contribution 

Z zero = exp /?y^ ero (-m + ^ )^™ > ( 26 ) 
and the Z' contribution is 

Z ' = nil ^(z^K + z^ + ^pfe + m)]" 1 . (27) 



3.1 The spin sector 



It is worth to notice that until this moment we have not used an explicit representation of 
the DKP algebra. In this point we specialize in the spin sector of DKP theory, thus we 
calculate every contribution in l|25() . We first compute the 1271) Z' partition function, to 
make it we use an explicit 5x5 representation (|58l) (sec Appendix) of the DKP algebra, 
thus we obtain 



In, 



z ' = -v E /tS ln {™ 3 [k + ¥? + - 2 ] } 



(2*) 



(28) 



where uj = \/p 2 + m 2 , and ignoring the /3-independent constant, we can rewrite 1)28(1 to 
show explicitly the existence of particles and anti-particles, thus 



d 3 p 



(2tt) 3 
d 3 p 



d 3 p 



In 



'2tt) 3 



(2tt) 3 

ln{/3 2 [oj 2 n + (oj + ^ 2 ]} . 



[p 2 [J* + (u- M ) 2 ] } + 



(29) 



To calculate the sum on n in the last two terms we use the standard results founded in 
the literature (EIGEjj thus 



\tlZ' 



-V 



2V\n((3)Y, 



'uj + In [1 

d 3 p 

(2tt) 3 ' 



In 1 



+ 



(30) 
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where the last term is an infinity /3-dependent quantity which has to be cancelled by the 
N((3) constant, thus we set it to be 

lnN(0)=-2Vln(f3)Y / Jj!^ (31) 
in total agreement with the literature when is studied the complex or charged scalar field 

0EDEH]. 

The zero mode contribution 

Now, we show the existence of the zero mode sector for the spin DKP field. The zero 
mode equation is obtained from the action given in (|21|l . 

(D + /^ )V = (32) 

by setting the fields ip to be independent of (r, x), thus 

( - m + /i/3°) i, zero = (33) 

then the determinant of the matrix y—m + /i/3°) must be set to zero for getting non-trivial 
solutions, 

det (-m + nJ3°) = -m 3 (m 2 - fi 2 ) = 0. (34) 

It is easy to see that we get non-trivial solutions when the chemical potential reaches the 
values ±m. The field V'zero is a five complex component column matrix, 

ipzero = (a, b, c 1 , c 2 , c 3 ) T . (35) 

Solving the matrix equation above, we find 

h fib mb 

c k = , a = — or a = — . (36) 
m fi 

Then, using the above relations we find that zero mode contribution of the partition function 
is for /i =/= ±m 

Z zero = -fiV(m 2 - M 2 )C 2 (37) 

with £ being an arbitrary real parameter. 

Finally, the full partition function for the spin sector of DKP field is 



InZ 



P V{m 2 -ii 2 )e-V f -0^ + m(l - e"^""^) + ln(l - («+/•))] (38) 



reproducing exactly the partition function of charged scalar field [Sj when we take the con- 
tribution of the zero mode into account. The momentum integral is convergent only if 
< m. The parameter £ that appears in the final expression for the partition function is 
not determined a priori, it should be treated as an variational parameter which is related to 
the charge carried by the condensed particles. Then, at fixed (3 and /i, In Z is an extremum 
with respect to variations of such a free parameter, 

^ = -2/3y(m 2 - M 2 K = (39) 

which implies that £ = 0, unless \fj,\ = m and in this case £ is undetermined by this variational 
condition. Thus, when < m the contribution of the zero mode to the partition function 
is cancelled. 
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3.1.1 The Bose-Einstein condensation 



By completeness, we briefly summarize JS] some characteristic of Bose-Einstein condensation 
in the case of a relativistic Bose gas. The question which we make is this: What is the 
requirement for BEC to take place at relativistic temperatures (i.e., T 3> to) and what is 
the nature of the phase transition. 

For \fi\ < to, the £ parameter is set to zero in (|38|) . and the charge density is given by 



1 /d\nZ\ f d 3 n ( 1 1 



(3V V dp J J (2tt) 3 yeP^-ri - 1 ^+^-1 



(40) 



as the chemical potential satisfies \p\ < to, thus, the charge densities of particles and anti- 
particles are non-negative. Note that (|40|l is really an implicit formula for p as a function 
of p and T. For T above the some critical temperature T c , one can always find a p such 
that l|4U|) holds. If the density p is held fixed and the temperature is lowered, p will increase 
until the point \p\ = to is reached, thus, in the region T > T c 3> to we obtain [9] 

\p\ « 1 toT 2 (41) 

When |//| = to and the temperature is lowered even further such T < T c , the charge 
density is written as 



1 fdluZ 



(3V V dp 



= Po + p*(j3,p = m), (42) 

/i— m 



where po = 2to£ 2 is a charge contribution from the condensate (the zero-momentum mode) 
and the p*((3,p = to) is the thermal particle excitations (finite-momentum modes) which is 
given by (|40|) with \p\ = to. 

The critical temperature T c , in which the Bose-Einstein condensation occurs, corre- 
sponds when \p\ — to is reached and determined implicitly by the equation 

p = p*(0c,V = m), (43) 

which implies that 

t c = (M)"\ m 



TO 

At temperatures T < T c , 1421) is an equation for charge density p — po of the p ^ states, 

P - Po = \ mT 2 , (45) 
so that the charge density in the ground state (p = 0) is 

(46) 





' T ' 












T 
.a . 





as the temperature is lowered the fraction of charge in p = state increase to unity. We 
note that l)44|l leads to the important result: the necessary condition for that any ideal 
Bose gas of mass to Bose-Einstein condense at a relativistic temperature (i. e., T c 3> m), 
is that p 3> to 3 . Additional literature about Bose-Einstein condensation can be found in 

[1I22EH1E1 
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3.2 The spin 1 sector 



Now we study spin 1 sector of DKP theory, thus we calculate every contribution in H25|) . 
We first compute the (|27H Z' partition function. For this purpose, we consider an explicit 
10 x 10 representation (|59"jl (see Appendix) of the DKP algebra, thus we get the following 
contribution 



In Z' = -Vj2j 7^3 ln |™ 4 [K + *M) 2 + 



(47) 



and more one time by ignoring the /3-independent constant, we can rewrite 147(1 to show 
explicitly the existence of particles and anti-particles, thus 



n 

9~ 2^ 



d 3 p 

d 3 p 



W 
~2~ 



(w-/i) 2 I 



In- 



(48) 



Here, we can see the six degree of freedom corresponding to two complex massive vector 
fields or Proca's complex fields. Carrying out the sum in both last terms and after some 
manipulations [511201, we nn d 



In Z' 



-W 



rf 3 T 



( .^[^ + In(l-e-^))+ln(l-e-«^)) 



(2tt)* 



(49) 



the last term is an infinity quantity /3-dependent, and, as in the spin sector, it will be 
cancelled by the N(f3) contribution, thus with all security we set it to be 



(50) 



The zero mode contribution 

In this case, we also show the existence of the zero mode sector for the spin 1 DKP field. 
The zero mode equation is obtained from the action given in (|21[1 . 

(d + /i^°)v = (51) 

by setting the fields tp to be independent of (r, x) , thus 

( - m + M/3°) i> zero = (52) 

then o determinant of the matrix (— m + must be set to zero for getting non-trivial 
solutions, 

dct (-m + /i/9°) = to 4 (to 2 -a* 2 ) 3 = . (53) 

It is easy to see that we get non-trivial solutions when the chemical potential fi reaches the 
values ±m, and in this representation the field ip zer o is a ten complex component column 
matrix, 

Vwo = (V, ^\i>\ V' 6 , i>\ ■ (54) 
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By solving the above matrix equation, we find 



i) T = t/> 4 = V 5 = ^ 6 = 

(55) 

! = /Z^ ^ 2 = ^ ^ 3 = Hip 9 

m m m 

By using the above relations we find that zero mode contribution of the partition function 
is for fi ^ ±m 

Z ze ro = ~P V(m 2 - /i 2 ) (fi + e 2 + it) (56) 

with £fc being an arbitrary real parameters related to the spatial components to the zero 
mode DKP field, i.e, ip 1 , ip 2 , t/j 3 . 

Finally, the full partition function for the spin 1 sector of DKP field is 

In Z = -(3 V(m 2 - y 2 ) + g + + (57) 

(3oj + ln(l - e-^-^ + ln(f - e ~^+^ . 

The parameters are related to the degree of freedom of the spin 1 DKP field, i.e, 
three degrees of freedom of the massive charged vector field. The zero mode sector is again 
absolutely necessary to reproduce the BEC as the spin zero sector; however, the critical 
temperature is the same for every degree of freedom. 

4 Conclusions and remarks 

The free massive DKP theory at Finite Temperature (FT) is equivalent to both complex 
massive scalar field and complex vector field theories at FT. It reproduces the relativistic 
Bose-Einstein condensation in both sectors, where we show in a clean and elegant way the 
zero mode existence and its contribution to BEC. 

The perspectives to follow are to study the DKP field coupled to the quantized electro- 
magnetic field and the implementation of the renormalization process at zero temperature 
which will allow to extend the analysis at Finite temperature. And, Bose-Einstein conden- 
sation in curved space-times |29j using the DKP theory minimally coupled. Advances in 
this directions will be reported elsewhere. 
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Appendix 

Spin and Spin 1 representation 

We use the following representation for the spin sector of DKP algebra 
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8° = 



















, p = 




-1 
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, P 3 = 
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(58) 



and for the spin 1 sector of DKP algebra 
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